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Abstract 

A popular experiment of the undergraduate physics laboratory is the study of a mass attached 
to a spring oscillating vertically. This experiment in its most general form, exhibits complex 
dynamics due to the simultaneous action of gravitational, elastic and torsional forces. In addition 
the friction of the air against the mass introduces further difficulties in the study of the dynamics. 
The objective of the article is to present a didactic method to study harmonic and non-harmonic 
motion in the spring-mass experiment in a comprehensive but still manageable level. The system 
is at first approximated as an harmonic oscillator. Significant deviations from the simple harmonic 
motion model are described and discussed. In particular this includes the observed parametric 
oscillation behavior. The study of the parametric behavior is done via the analysis of the motion 
waveforms. 
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I. INTRODUCTION 



An experiment featuring a mass attached to a spring is a popular experiment worldwide. 
In many instances it is the basic experiment used to introduce first year physics students 
to the study of harmonic oscillations. This indeed was the case in our first year physics 
undergraduate course 1 . A key advantage of this experiment is that the apparatus is simple 
and relatively inexpensive. A downside, however, is the number of papers indicating that 
this experiment is not simple at all^~-. 

Gravitational (pendulum), elastic (spring) and torsional forces^ 1 ^ generates numer- 
ous and complex physical phenomena in the motion of the spring-mass system. Examples 
are: multimode operation 8 , parametric instability and energy transfer between the spring- 
bouncing mode and the pendular-swinging mode. A real spring-mass system behaves as a 
simple harmonic oscillator-as treated in textbooks-only in the specific conditions: (i) the 
mass of the spring is negligible with respect to the attached mass, (ii) the system oscil- 
lates far from resonance, that is the ratio between the spring Uk and pendulum u p mode 
frequencies is far from 2, (iii) the system featuring a spring with an exact axial stretching. 
In general the spring-mass behavior in a laboratory environment opened to many possible 
choices turns out very complex to study. In fact, it behaves as a parametric oscillator rather 
than an harmonic oscillator. This is due to the fact that the spring and pendulum motions 
are non-linearly coupled. In this case the oscillation becomes unstable when the frequencies 
of the two modes are resonant^, j. e . u) k /u) p = 2. 

We believe that a comprehensive experimental treatment of the general spring-mass sys- 
tem is needed. This is mainly due to the nature of the laboratory environment where students 
tend to operate the system in a variety of possible configurations. The main punch line of 
the paper is to show that the study can be analyzed via the motion waveforms obtained 
with a set of different masses. Different masses induce a variation of both the spring and the 
pendulum angular frequency (via the pendulum length variation) and thus of the frequency 
ratio 



where k is the spring constant, m is the appended mass, £q is the rest spring length and g 
is the gravity acceleration constant. A suitable change of the masses allow to span all the 
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behaviors from below resonance to beyond resonance. 

In running the experiment, students have to face the gap between theory and practice, 
have to look at a complex multi-effect motion, have to master experimental techniques, 
have to treat data statistically and, possibly, have to purge the experiment from the the 
many "nuisances" (with the help of a guide). They are forced to develop the ability to 
manage unexpected experimental observations. The complexity of a real behavior requires 
a strategy for singling out the motion components, relaxing their interrelations so to treat 
separately the phenomena and, afterwards, to set the complete picture. In short, the work 
on this experiment turns out to be a research (guided) activity. The physics content of a 
parametric behavior must be approached and the basics must be enlighten. 

The theoretical work starts from the textbook idealized motion equation 

d 2 z(t) C dz(t) k , , . „ . , . 

at z m at m 

The experimental goal is to test the solution 

z(t) = Zee -1 * sin{u)t) (3) 

where 

^ 2 = --(^) 2 = ^-7 2 . (4) 
m \2mJ 

In these equations, 7 = C /2m is the damping constant and C is the friction coefficient. The 
damping time is r = I/7. In the above equations, the gravitational force mg is opposed to 
the equilibrium spring force k£o (hence —k£o + mg = 0) and sets the rest position z(0). The 
mass of the spring is assumed zero, while it is not at all. Friction is assumed proportional 
to the velocity, while air resistance depends on the object shape and, in fact, it is not linear 
with velocity. The investigation proceeds trough the following steps: (a) test of the Hooke's 
\F\ = \k Az\ law; (b) test of u> 2 = kjm law; (c) measure of the damping time r and (d) test 
of the sinusoidal solution of Eq. ([3]) . 

Here we present the study of the free spring-mass oscillations. A complete investigation 
of the spring-mass system requires also the study of the forced oscillations^. This study is 
meant as a deeper grounding for those conducting a widely used experiment and a possible 
topic for undergraduate students in open ended investigations or small individual research 
projects. 
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II. THE SETUP 




FIG. 1: Layout scheme: a meter and light projector are used for measuring the spring extension, the 
sonars detect the motion, the digital data acquisition tool LPRO - a Vernier Loggerpro- acquires 
and sends the coordinates z(t), x(t) to the computer. The top arm is bent for avoiding back 
reflection of the sonar signal. 

The apparatus has been kept as simple as possible, just a spring and a weight oscillating 
vertically as shown in Fig. [TJ The pivot is done by a short piece of wire. Sonarl and sonar2 
detect z(t) and x(t) coordinates. The system is surrounded by plastic foam to absorb the 
direct sonar waves not reflected by the mass-blob. This is a pile of 20 g metal disks on a 
support, to which we added as much plasticine as necessary to obtain precise values of total 
weight. The mass moves around in a almost random way owing to a mass wobble motion and 
the spring transverse oscillations (with wide pendular and/or vertical oscillations). Because 
of these spurious motions, at the base, the support holds a diskette of 70 mm diameter to 
reflect efficiently enough the sound waves sent by the bottom sonar. The diskette diameter 
value is a compromise between a larger value required for ideal sonar detection and the 
smallest possible one required to minimize wobble motion. Spurious motions are kept under 
control by carefully stretching the spring along its vertical axis or, alternatively, shifting 
the mass sideways. The resulting waveforms are therefore as simple as possible. A generic 
system motion is a composition of the two spring and pendulum oscillations plus the mass 
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FIG. 2: Scheme of the sound waves paths to and fro sonars-masses. Sonars are suited to follow 
the in front movements. A side movement is detected as a length difference between the initial 
and actual sonar-mass distance, the side movement is seen as a kind of projection of the initial 
sonar-object line along the actual sonar-object line. Sonarl detection is indifferent to the plane of 
pendulum oscillation, the inverse happens for sonar2. 

wobbling. The bottom and side sonars detect the composite motion, see Fig. |2j The bottom 
sonar sees the transverse (pendulum) motion with a very reduced amplitude (its projection 
onto the z axis) and the converse happens for the second sonar. The distortion introduced 
into one motion by the other motion was tested by forcing the system into a pure motion. 
Just for this test, the pendulum motion was obtained appending the mass at a wire, the 
vertical motion was forced inserting a thin metal stick vertically inside the spring. 

We will see below that the pendulum oscillation plane is not stable. This fact does not 
affect the detection by sonarl, whereas it makes sonar2 useless when the angle variation 
is large, but still helpful when the angle variation is small. The motion analysis is based 
mainly on data from sonarl. 

Typical mass and spring constant values were m spr i ng ~ 4.8 g, k ~ 8 N/m. The following 
effective masses were appended: 40, 45, 50, 55, 60, 70, 80 g. The effective mass is the sum 
of the appended mass and l/3rd of the spring mass^&ii, as explained below. In fact, this 
experiment also requires the consideration of the non-negligible mass of the spring. 

III. TEST OF THE HOOKE'S AND uj 2 = k/m LAWS 

Hooke's law \F\ = \k Az\ test was performed by appending a set of different masses and 
measuring the spring elongations with a meter. The spring constant value was found with 
a precision of 0.1 percent. 
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The u 2 = k/m law test was done appending a set of different masses and measuring the 
period T by a chronometer. Masses range from a top value limited by the spring damage 
and a bottom value determined by excessive vibrations. The goal was to verify the above 
relation transformed into 



k 



m 



4tt 2 



• T 



(5) 



The graph of the mass versus T 2 showed the expected straight line, but with a negative 
intercept at the y-axis, see Fig. [3j This test posed "dramatically" the problem of the 
nonnegligible mass of the spring. The class discussion about the physics lawi led to conclude 
that the spring mass had to be considered and the equation was rewritten as 




FIG. 3: Graph of eq. [5] showing the negative intercept at the y axis. 



m = ■ T 2 - m . (6) 

The intercept gave the expected value of a third of the spring mass to be added to the 
dynamic system^- 1 m = (1/3) m spr i ng . In doing the above test, students noted that the 
lighter the masses the more irregular the oscillation motion and the more scattered the 
results. They realized that the system did not show the expected simple harmonic oscillation 
and posed questions. The new physics of the parametric oscillation made its first apparition 
in our undergraduate Lab. 

An estimate of the decay time r is done measuring with the chronometer the time at 1/e 
oscillation amplitude reduction. Different masses produced different r as expected from the 
equation r = 2m/ C. 
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IV. MEASUREMENTS ON HARMONIC AND PARAMETRIC MOTION 



After the above tests we passed to the motion analysis using the motion waveforms. 
Students had to test motion Eq. (j2j) via the analysis of the detected waveform fixed at the 
computer screen. They operated with different object masses following their own project. 
Typical waveforms obtained with the set of 70 g, 60 g, 50 g and 40 g masses are shown in 
Fig. HI Waveforms show that the spring-mass system has a substantial vertical harmonic os- 





FIG. 4: Waveforms from top to bottom refer to 40 g, 50 g, 60 g and 70 g masses. The system 
passes from near resonance, with 40 g mass ujk/ujp ~ 2.04, to out of resonance, with 70 g mass 
Wfc/wp ~ 1.67. The waveforms evolve from complete modulation to classical damped sine-wave. 
The comb-like crests along the signal envelope, evident in the first signal, are generated by the 
peculiar detection of the transverse motion by sonarl. 

cillation mode with relatively heavy masses, but two modes with light masses. The presence 
of the two modes modulate the envelope. The direct observation of the motion shows that 
in certain configurations (i.e. proper ratio between the spring-bouncing and the pendular 
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frequencies) the vertical motion disrupts in few cycles. The plane of the pendular oscillation 
changes in an apparently random manner. 

The harmonic and non-harmonic behaviors refer to different physical phenomena. Fig. 
clearly shows the coupling between the two spring and pendulum oscillations. It shows also 




FIG. 5: Sketch of the applied force with coupling between spring-bouncing and pendular-swinging 
motions. Note that the mass does not align with the spring axis. 

that resonance between the two motions occurs when the ratio between two motion periods 
Tfc (spring) and T p (pendulum) is 1/2. In the pendular-swinging motion the stretching 
force applied by the mass to the spring causes a double oscillation in a single period of the 
pendular motion. This force oscillation makes the "rest spring length" to oscillate about its 
initial rest value £q 

e (t) = e + se(t). (7) 

The coupling between the pendular and spring oscillations induces a parametric instability 
between the two modes^*^ when the resonance condition is met. The degree of instability 
and the mode energy transfer depend on the angular frequency ratio of Eq. (TjQ). This 
is easy to check by exciting the pendulum motion with a lateral displacement. One sees 
that vertical spring oscillations are generated at the expense of the pendulum motion. The 
converse happens with an initial vertical stretching. Fig. [6] shows the induced vertical signal. 



are given as^ 



In literature, equations describing nonlinear coupled pendulum and spring oscillations 

^-^- + u 2 p x(t)=cx(t)z(t) (8) 
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FIG. 6: This waveform of the vertical oscillation comes from the energy transfer from a pendulum 
oscillation carefully started by a lateral initial displacement. The "fuzzy" shape of the waveform 
depends on the simultaneous presence of vertical and transverse oscillations, as seen by sonar 1. A 
70 g mass was used in this test. 

d 2 z(t) 9 . . x 2 (t) 

where the coupling constant c is a function of the system parameters 2 . The way to obtain 
these equation is through the Lagrangian formulation in which friction force terms are not 
considered. This system can be solved in a particular case analytically 2 and in the gen- 
eral case numerically^. The onset of the parametric instability with a complete conversion 
between the two oscillations occurs when the resonance condition is met. 

Modulated waveforms, in which there is only partial energy transfer between modes, were 
observed up to values of the ratio Uk/co p around 1.8. We suggest that the presence of the 
disk with its wobble oscillation might produce a widening of the resonance condition for the 
mode exchange. These observations are supported by the analysis of the frequency spectrum 
given by the waveform FFT (Fast Fourier Transform), shown in the frames of Fig. [7J In 
each frame, the largely dominant line is that of the oscillation. Other lines correspond 
to pendular motion; due to the peculiar detection of the transverse motion by sonarl (see 
Fig. [2]) the relative second harmonic frequencies are particularly evident. The spring and 
the pendulum frequencies resulted shifted of a couple of centimes of Hz, as foreseen by 
theory- 2 *. We checked that the very small lines at high frequencies are due to the disk wobble 
vibrations. 

Incidentally, students had to measure the frequency from the waveforms as a basic test 



of Eq. (j4]). In that connection the concept of frequency spectrum and the use of a FFT tool, 
giving immediately the frequency, was introduced. 
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FIG. 7: FFT signals derived from 40g, 50g, 60g and 70g waveforms. FR is the actual value of 
the frequency ratio. The frequencies of some evident spectrum lines are indicated. The neater the 
waveform the neater the FFT. 

Modulated waveforms were observed with relatively heavy masses in the case of mixed 
(vertical and transverse) excitation as, for example, that presented in Fig. |HJ The setting of 
an initial x(0) ^ causes the onset of the exchange between the two modes, as predicted by 
the coupling term in the right hand side of Eqs. OH]) and In this case, after few oscillation 
cycles, the system sets itself in the peculiar stable motion shown in Fig. [9j described also in 
re£2. The FFT spectrum shows many relatively intense lines, as in the second frame of Fig. [Sj 
They are generated parametrically by the nonlinear coupling between the two fundamental 
oscillation modes. 

General considerations on the parametric observed behavior are reported below. 
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FIG. 8: This waveform was obtained with 80g mass, ratio uJk/^p — 1-63, and small diagonal 
displacement. The exchange between the two modes was visibly clear. The frequency spectrum, 
second frame, shows many lines besides the spring frequency at 1.55 Hz and the very small pendular 
frequency at 0.95 Hz. Two of these lines correspond to second harmonics of pendulum and spring 
frequency, respectively. Other evident lines at higher frequency correspond to wobble motions. 




FIG. 9: Motion figure normally assumed by masses heavier than 60g. 

• Both strong vertical and lateral mass displacements lead to substantially similar results 
in terms of waveforms and spectra. This is not the case with small displacement: the 
lateral one moves both the pendular and the vertical force, the vertical displacement 
does not switch the pendulum force. 

• The mixed excitation produces complete energy transfer between the two modes also 
at Uk = 1.8 u p . The spectra are neater than in pure vertical excitation. 

• The pendulum oscillation plane changes at each mode exchange. 
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A. The connection with the optics 



About the parametric behavior of our system, it is worthwhile reporting that the para- 
metric instability is extensively exploited in optics to produce and amplify laser waves at any 
requested frequency^ 3 . A pump wave u\ launched across a non-linear crystal generates a 
signal wave Co>2 together with an idler wave Ui. The condition that must be satisfied is 

uj x = uj 2 + ooi (10) 

In our spring-mass system the spectrum shows clearly the so-called idler frequency Ui 

Ui=U k -U p . (11) 



In many tests, this oscillation appears very small so it must be searched in the spectrum at 
a known position. As an example, the idler frequency line obtained with a mixed excitations 
is shown in Fig. [TUJ 
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FIG. 10: The idler frequency line at 0.63 = (1.64 — 1.01) Hz from the FFT graph of the waveform 
generated by vertical excitation with an effective mass of 70 g and recorded by sonar 2. The same 
result was obtained with the transverse (pendular) excitation. 



V. MEASUREMENT OF THE DECAY TIME r FROM WAVEFORMS 

Following the program on the study of harmonic oscillation, the waveforms far enough 
from resonance were analyzed. Students obtained a value for r from the displayed waveform 
by measuring the time from the start to 1/e of the initial amplitude. The decay time 
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resulted much longer than the value previously quoted by direct measurement. Moreover, the 
oscillation persisted for a very long time at the tail. The waveform should fit the theoretical 
solution of Eq. (J2J) which is the mathematical function f(t) = a(0) ■ exp(—t / r) ■ sin(u + (p) + b 
by substituting a(0),r, v obtained in the initial direct measurement. The fit turns out to 
be impossible. Reasonably good fits between the theoretical curve and the waveform are 
obtained within sections of the decaying waveform. As an example we report the results 
obtained with 70 g mass. The fit using the first 10 s curve gives r ~ 12 s, the fit in the 
range At = — 20 s gives r ~ 17 s, the fit in the range At = 10 — 30 s gives r ~ 20 s, 
the FIT in the range At = 20 — 60 s gives r ~ 30 s. The oscillation frequencies varied in 
the different curve sections in agreement with the r variation: longer r , higher frequency. 
Moreover, the stronger the excitation, the shorter the measured r in the first curve section. 
The persistence of the oscillation exceeded 100 s. The fitting curves showed a tail amplitude 
thinner than the recorded waveform. 

To account for these observations, we will consider the energy exchange between the 
vertical and horizontal modes. The decay of the oscillation waveform detected by the bottom 
sonar has two ingredients: one is the energy transfer to the pendulum mode, the second is 
the lower friction at lower oscillation amplitude. The initial fast decay might reasonably 
be due to the energy transfer, while the slow decay at the tail is simply due to the strong 
reduction of the air friction with a slow-moving disk (at that point there is no exchange 
between the two oscillations). 

The model for friction 7a; = (C/m)x states that the decay time r is proportional to 
the appended mass. Measurements result in agreement with that assumption at an extent 
better than 80 percent^. 

VI. CONCLUSION 

The dynamics of a generic spring-mass system moving vertically is equivalent to that 
of a parametric oscillator with energy exchange between the spring-bouncing and pendular- 
swinging oscillation. The spring-mass physics system used in a first year laboratory course for 
the investigation of harmonic oscillations has to be carefully designed in order to reproduce 
the characteristics of the idealized system treated in textbooks. In particular, the system 
has to be far from resonance and must feature an exact axial spring stretching. 
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The study of experiments constrained along a narrow pre-defined path has a poor didactic 
content. It limits the thinking of the student as a scientist: developing a model and testing 
it. Instead, for a student to grow, it is paramount to investigate physical phenomena in open 
environments. We believe that the spring-mass experiment should be approached in class 
as an harmonic oscillator. However the laboratory environment should be left unstructured 
granting freedom to the students to build their own system by choosing spring types and 
different masses. This approach will be an effective teaching tool exposing the students to 
the dense physics contents of the parametric behavior of the spring-mass system. 

The investigation of the motion of the spring-mass experiment via the analysis of the 
motion waveforms has the potential to enlighten many physics phenomena present in that 
system: motion instability, mode transfers, resonance, the passage from parametric to har- 
monic oscillation, the friction variation with oscillation amplitude. 

The main contribution of this paper is to show that the spring-mass experiment is man- 
ageable and well within the reach of the curriculum of a first year physics course using the 
waveform detection and its analysis. The complexity and the gap between the ideal theoret- 
ical case and the real case makes this experiment a very useful tool to stimulate the thinking 
process of an aspiring physicist. 
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